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I.  INTRODUCTION 


Cable  shields  and  cable  conduits  have  been  widely  used  as  part  of  the 
integrated  electromagnetic  pulse  (EMP)  hardening  design  for  aircraft,  missiles, 
and  ground-based  systems  (Ref.  1).  When  cable  shields  are  used  over  a  long 
distance,  they  are  generally  periodically  grounded  to  their  immediate  outer 
shield  through  bonding  straps,  clamps,  screws,  etc.  The  periodic  grounding 
provides  mechanical  rigidity  of  the  shields,  reduces  the  electrostatic  hazards, 
and  at  the  same  time  might  increase  the  shielding  effectiveness  against  pene¬ 
tration  of  long-wavelength  disturbances.  However,  for  the  broad-band  EMP,  the 
question  arises  whether  the  periodic  grounding  will  improve  or  degrade  the 
overall  shielding  effectiveness.  In  this  report,  this  question  will  be 
answered  from  a  rigorous  theoretical  analysis.  The  analytical  results 
will  help  the  system  engineers  to  determine  the  optimum  grounding  arrangement. 

Cable  shields  generally  are  not  perfect.  EMP  can  penetrate  either 
locally  (such  as  through  an  aperture,  a  connector,  or  one  end  of  the  shield, 
see  Fig.  1)  or  distributedly  (such  as  through  diffusion  or  uniformly 
distributed  apertures,  see  Fig.  2).  These  two  distinct  cases  will  be 
separately  discussed:  discrete  excitations  in  Section  II  and  distributed 
excitations  in  Section  III.  Finally,  the  results  obtained  in  Sections  II 
and  III  will  be  summarized  in  Section  IV,  and  several  practical  engineering 
examples  as  to  how  the  results  can  be  used  will  be  presented. 
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II.  DISCRETE  EXCITATIONS 


In  this  section,  the  outer  shield  of  a  double-shield  configuration 

on  Figure  1  is  assumed  to  be  perfect,  except  at  some  isolated  locations 

where  the  EMP  can  penetrate.  The  current  and  charge  density  distributions 

on  the  transmission  line  formed  by  the  two  shields  due  to  the  penetration 

from  a  localized  location  can  be  calculated  by  solving  the  problem  depicted 

in  Figure  3.  The  distributions  due  to  the  voltage  and  current  sources  of 

the  localized  penetration  can  be  separately  considered.  The  two  problems 

of  Figure  3a  can  be  decomposed  further  (  Fig.  4  )  so  that  only  the  problem 

shown  in  Figure  3b  (which  is  redrawn  as  Fig.  5  with  Vs  *  V?,V*  or  Vs  and 

I  L  O 

s  s  s  s 

I  ■  I^,^  or  I0)  need  be  solved.  The  coupling  to  the  wires  inside  the 
inner  shield  can  then  be  calculated  by  multiplying  the  appropriate  transfer 
impedance  or  charge  transfer  frequency  of  the  inner  shield  by  the  current 
or  charge  density  distributions.  The  time  variation  of  exp(jiot)  will  be 
assumed  and  suppressed  throughout  this  report. 

At  sufficiently  low  frequencies  the  TEM  mode  is  dominant  and  the 
voltage  (V)  and  current  (I)  distributions  along  the  line  shown  in  Figure  5 
can  be  calculated  by  solving  the  following  pair  of  equations: 

S  +  K1  • 0  <l) 

V  *  0,  for  z  >  z  (2) 

o 


dl 

dz 


Y*  +  Y 
1  *d 


l  <5(z  -  nd) 


where  and  Yj!  are,  respectively,  the  series  impedance  and  shunt  admittance 
per  unit  length  of  the  line  and  is  the  admittance  of  each  strap. 

From  Floquet's  theorem,  the  solution  of  the  periodic  Equations  1  and  2 
has  the  following  form: 

00 

V(z)  ■  £  V  exp[-j(k  +  2nr/d)z]  (3) 

n«-°°  n 

00 

I(z)  •  l  Inexp[-j(k  +  2nw/d)z]  (4) 

n*-® 


12 


Figure  5.  The  resultant  tra 


Substituting  Equations  3  and  4  into  Equations  1  and  2  one  has 


ln  “  j  Vn(k  +  2n1T/d)/z[ 


(5) 


where 


V 

n 


-  Y  Z' 
d  1 


YjZ'  +  (k  +  2n7r/d) 


2  d 


00 

v°  =  y  v 

*•  n 
n=-°° 


(6) 


(7) 


By  combining  Equations  6  and  7,  one  finds  that  the  dispersion  relation 
between  k  and  /-Y|Z|  is  (also  see  Ref.  2) 


cos  kd  =  cos(d/-Y^Z j) 


ZIYdd 

2d/-YjZ| 


sin(d^-Y|Z|) 


(8) 


In  Figure  6,  cos  kd  is  plotted  versus  d/-Y jZ|  (=  ood/c)  for  various 

qg  (*  Z|Y^d) .  (Generally,  Z| ,  Y|  and  Y^  are  complex  values.  To  have 

real  and  positive  d/-Y' Z '  and  Z'Y.d,  one  can  assume  Zl  =  Z'  ,  Y'  =  Y' 

1  1  Id  1  cl  1  cl 

and  Y^  a  (juL^)  ,  where  Z^  and  Y^  are  purely  imaginary  and  are 

respectively  the  values  of  Z^  and  Y^  of  a  perfect  double  shield,  Ld  is 

the  inductance  of  each  bonding  strap).  In  the  figure,  one  clearly  sees 

the  passband  and  stopband  structures.  The  stopbands  are  generally  broader 

for  a  larger  qg,  especially  at  the  lower  wd/c  region.  This  is  reasonable 

because  a  larger  qg  means  a  better  grounding.  At  the  higher  wd/c  region, 

the  stopbands  become  narrower  because  the  high-frequency  disturbance  does 

not  see  the  presence  of  the  bonding  straps.  Also,  to  have  the  proper 

propagation  and  decaying  constants  for  z  >  z  ,  one  should  restrict  the 

o 

k-value  to  have  a  positive  real  part  and  a  negative  imaginary  part.  The 
imaginary  part  of  k  determines  how  fast  the  disturbance  decays  when  it 
propagates  away  from  the  penetration  point.  A  plot  of  |lm(kd)|  versus 
frequency  for  various  qg  is  given  in  Figure  7.  Figure  7 

obviously  shows  that  a  larger  qg  gives  a  larger  [ Im(kd) |  (which  means 
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Figure 


Figure  7.  The  decaying  constants  of  bonded  transmission  lines  versus  frequency  for 


that  the  disturbance  decays  faster)  and  thus,  is  more  desirable.  Actually, 
at  the  low-frequency  limit  where  d  =  wd/c  <<  1,  one  can  easily  show 
that 


tu 


i 


L_1(l  +  qs/2) 

(9) 

rX[l  +  L[d/(2Ld)] 

(10) 

where  is  the  series  inductance  per  unit  length  of  the  double  shield. 
Using  Equations  5  and  6  in  Equations  3  and  4,  one  has 


V(z)  -  l 


(Y|Z|  +  k‘) 


0  n=-»  Y|Z|  +  (k  +  2nir/d) 2 


-j  (k  +  2nir/d)  z 


YdZld 


cos(kd/2)cos  £  (  (m  +  l/2)d  -  z)  /-Y|Z^  j 
cos(d/-Y|Z^/2) 

sin(kd/2)sin[(  (m  +  l/2)d  -  z)/^Z|  ] 


sin(d/-Y^Z|/2) 


(11) 


for  md  <  z  <  (m  +  l)d 


I(z)  *  Vn  l 


j(YjZ^  +  k  )  k  +  2nir/d  -j  (k  +  2nir/d)z 


0  n=-“  YjZ|  +  (k  +  2nff/d) 2  Z1 


e-jkd(m  +1/2)  x 


-j  VQd/-Y|Z|(k2+YjZj)d 

<dVi)zid 

|Sin(kd/2)cos  (m  +  l/2)d  -  z)  /-Y|Z|  ] 


sin(d/=Yp^/2) 
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cos(kd/2)sin  F  ( (m  +  l/2)d-  z)/-Y|Z'  1 
+  j  - - - -  - - —  J 


COS  (d/^YjZj/2) 


(12) 


for  md  <  z  <  (m+  l)d 


To  obtain  Equations  11  and  12,  one  has  used  the  dispersion  relation  8  and 
several  series  summation  formulas  in  Reference  3.  With  the  expression  for 
I(z),  one  can  calculate  the  total  charge  per  unit  length  of  the  double  shield 
from 

Q’(z)  =  -  -r-  f-  I(z)  (13) 


As  is  obvious  from  Equations  11  and  12,  Vq  is  yet  to  be  determined 
from  the  boundary  condition  at  z  =  zq.  In  what  follows,  the  cases  with  the 
voltage  and  current  sources  will  be  considered  separately. 

1.  VOLTAGE  SOURCE 

s 

Applying  the  condition  that  V  =  V  at  z-zq,  one  immediately  obtains 
from  Equation  11  that 


V 

o 


~ZiYd  V 


(k  +Y^Z[)d 


eJkd/2  x 


cos (kd/2)cos [ (z  -  d/ 2) /-Y’ Z' ] 
o  11 


cos(d/-YjZ[/2) 


sin(kd/2)sin[(zo-d/2)/-Y|zJ 

sin(d/-Y|Z^/2) 


-1 


(14) 
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The  voltage,  current  and  charge  distributions  of  the  double  shield  can  be 

fully  expressed  by  Equations  8,  11,  12,  13  and  14  provided  Vs  is  known. 

For  the  problem  of  Figure  3b,  Vs  is  simply  V®.  However,  to  obtain  Vs 

(Note,  one  uses  Vs  *  vf  for  z  >  z  and  Vs  *  vf  for  z  <  z  )  for  the  problem 

1  o  zo 

of  Figure  4a,  one  needs  to  go  through  some  complicated  algebraic  manipula¬ 
tions  involving  Equations  12  and  14  and  the  relationships  given  in  Figure  4a. 
For  the  special  cases  of  zq  *  0  and  d/2  which  are  the  cases  to  be  considered 
in  the  following  discussions,  Vs  -  *  V®  *  V^/2.  can,  generally,  be 

obtained  from  the  short-circuit  current  I  at  the  outermost  surface  of  the 

sc 

double  shield  and  the  transfer  impedance  Z,^  of  the  outer  shield  via 


1 

sc  T1 


(15) 


Since  I(z)  and  Q'(z)  are  now  known,  one  can  begin  to  discuss  the 

coupling  to  the  wires  inside  the  inner  shield.  The  coupling  can  be 

completely  described  by  V'(z)  and  I'(z)  which  are,  respectively,  the 

s  s 

voltage  and  current  source  terms  of  the  transmission-line  equations  of 
the  inside  cables.  If  the  inner  shield  has  a  shield  transfer  impedance 
per  unit  length  and  a  charge  transfer  frequency  then,  and  1^ 

can  be  calculated  from  (see  Ref.  4): 


V^(z)  =  Z^Kz)  (16) 

r(z)  -  nT2Q’(z)  (17) 


With  the  above  expressions  one  can  then  define  an  effective  transfer 
impedance  per  unit  length  Z,^  and  an  effective  charge  transfer  frequency 
per  unit  length  for  the  double  shield  as  follows: 


-  v;(z)/isc 

(18) 

!TV(z)  -  KW/1BC 

(19) 
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Here,  the  subscript  "V"  is  used  to  indicate  that  the  quantities  are  for  a 
localized  voltage  source  which  is  generated  by  I  . 

The  two  quantities  and  are  z-dependent.  When  the  frequencies 
are  in  the  stopbands,  they  become  extremely  small  for  z  >>  zq  (see  the 
exponentially  decaying  term  in  Equation  12).  On  the  other  hand,  when  the 
frequencies  are  inside  the  passbands,  they  are  oscillatory  functions  of  z 
and  are  modulated  by  the  periodicity  of  the  bondings.  It  thus  appears 
that  in  order  to  better  quantify  the  coupling  to  the  wires  inside  the 
inner  shield,  one  should  apply  some  average  schemes  to  and  fi'  over 
Che  period  d.  A  natural  scheme  is  to  define 


where  nir  £  kd  <  (n  +  l)ir  when  mr  <_  d  <  (n  +  l)ir,  n  =  0,1,2,  ....  These 

two  average  quantities  are  calculated  to  be 


j(  ZT1 ZT2 

/  vv^l75! 


nT2ZTlYl 


/_Y '  Z ' 
11 


kd 


d 


-(Z]Y  d)(kd)exp(jkd/2) 

Tzm m  2  2~  x 

('CziYId)(k  +YizPd 

cos(kd/2)cos[(z  -  d/2)/-z;YH 
_ o _ 1  1 

cos(d/-ZjYj/2) 


sin(kd/2)sin[(zo  -  d/2) *0^1 
sin(d/-ZjY|/2) 


(21) 


where  n 

P 


1  for  Figure  3b  and  n 


P 


2  for  Figure  3a. 
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Z’  and  £2'  still  depend  on  where  z  =  z  is.  However,  no  matter  what 
TV  TV  o  _  _ 

zq  is,  one  can  easily  observe  that  the  normalized  |z^v|  and  (normalized 

with  respect  to  their  values  when  there  is  no  bonding,  which  are,  respectively, 

ZtiV^p^I)^  ZTl^T2Y]/np^  are  2reater  than  unity  for  some  regions  of 

/-'z-Y;  d  (=  wd/c)  and  smaller  than  unity  for  the  others.  That  is,  the 

bondings  improve  the  shielding  effectiveness  of  the  double  shield  at  certain 

frequencies  but  degrade  it  at  the  others.  Plots  of  the  normalized  |Z^,V| 

and  I £2'  , I  for  z  =0  and  d/2  as  functions  of  ^-zIYT  d  and  qr  (=  Z-fY.d)  are 
TV  o  x  x  p  1  a 

given  in  Figures  8  through  11.  In  the  figures,  the  curves  are  not  given 
for  the  stopbands.  In  the  stopbands,  the  normalized  quantities  are 
exponentially  decayed  away  from  the  penetration  point.  Since  the  ranges 
of  the  stopbands  increase  with  q^,  one  should  try  to  have  q^  as  large  as 
possible  for  better  shielding  effectiveness.  This  can  be  achieved  by 
making  Z^Y^  (=  L^/L^)  large.  Note  that  making  d  large,  although 
increasing  q_,  widens  o)d/c  and  hence  narrows  the  stopbands.  From  the 
figures,  one  also  sees  that  arranging  the  bondings  with  d  *  cnn/u  (where 
resonances  occur)  seriously  degrades  the  shielding  effectiveness,  and  thus 
should  be  avoided  for  the  important  parts  of  the  EMP  spectrum. 

2.  CURRENT  SOURCE 

s 

Applying  the  boundary  condition  that  I  =  I  at  z  =  zQ,  one  obtains 
from  Equation  12: 


V 

o 


j  I  (dY^Z^)Z^d 
dk-YJZJ  Ck2  +  Y’Zj>d2 


eJkd'2  * 


sin(kd/2)cos[(zQ  -  d/2) /-Y|Z^] 
sin(d^-Y[Z^/2) 

cos(kd/2)sin[(zQ  -  d/2)/-YjZ^] 
cos(d/-Y^Z[/2) 


(22) 
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K-d-*i 
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Figure  9.  Normalized  [ft'  J  when  z  =0  versus  frequencies  for  various 


Normalized  |  when  z  =d/2  versus  frequencies  for  various 


The  voltage,  current,  and  charge  distributions  of  the  double  shield  due  to 

the  current  source  shown  in  Figure  5  can  then  be  completely  described  by 

Equations  8,  11,  12,  13,  and  22,  provided  that  Is  is  known.  Here,  similar 

s  s 

to  the  voltage-source  case,  for  the  problem  of  Figure  3b,  I  is  I  • 

g  g  g  S  S 

To  obtain  I  (1*1^  for  z  >  zq,  I  «I  for  z  <  z^)  for  the  problem  of 

Figure  3a  or  4b,  one  has  to  solve  Equations  11,  12  and  the  equations  given 

s  s  s  s 

in  Figure  4b.  For  the  special  cases  of  zq  *  0  and  d/2,  I  ■  ^  “  IQ/2. 

From  these  current  and  charge  distributions,  one  can  calculate  the  voltage- 
and  current-source  terms  of  the  transmission-line  equations  of  the  inside 
cables  using  Equations  16  and  17. 

g 

The  current  source  of  the  shield  I  can  generally  be  estimated  from 

the  short-circuit  charge  density  Q'  on  the  outermost  surface  of  the  double 

SC  L 

shield  and  the  localized  charge  transfer  frequency  (which  has  a  dimension 
of  frequency  x length,  different  from  that  of  of  the  outer  shield  via 

(for  example,  see  Equation  38  of  Reference  5) 


Is  -  q'  niL 
o  vsc  T1 

(23) 

One  can  also  define  an  effective  transfer  impedance  Z ^  and 
charge  transfer  frequency  for  the  double  shield  as  follows: 

an  effective 

ZTI(z)  " 

(24) 

an(z)  ■ 

(25) 

where  the  subscript  "I"  is  used  to  indicate  that  the  quantities  are  for 
a  localized  "current"  source  which  is  generated  by  Q’^. 

Similar  to  Z ^  and  of  the  voltage-source  case,  Z,^  and  are 
also  z-dependent.  They  become  very  small  far  away  from  zq  in  the  stopbands, 
are  oscillatory  and  modulated  by  the  periodicity  of  the  straps  in  the 
passbands.  The  averaged  quantities  and  8,^  are  given  by 
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ZIYI  d 


“tA2'CJ;Y:  d 


(Z ' Y  ,d)  (kd) exp ( j  kd/ 2) 


-ZjY|  d  (k  +YjZ|)d 


sin(kd/2)cos [ (z_  -  d/2)/-Z’Y' ] 


sin(d/-ZjY|/2) 


cos(kd/2)sin[ (2  -  d/2) /-Z'Y' ] 


cos(d/-ZjY^/2) 


ZTI  and  also  depend  on  where  z = Zq  is.  However,  no  matter  where 


zq  is,  one  still  sees  that  Che  normalized  |ZTI| 


and  ft_T  (normalized  with 


respect  to  their  values  when  there  is  no  bonding,  which  are,  respectively, 
Q^lZ^/Cjiuiip)  and  >  are  8reater  than  unitY  for  some 

/-Z|Y[  d  (where  the  bondings  degrade  the  shielding  effectiveness)  and 
smaller  than  unity  for  the  others.  Examples  of  the  normalized  jz^l  and 
|^TII  for  zo=*d/2  as  functions  of  /-Z^Y^  d  and  qg  are  plotted  in  Figures  12 
and  13.  Similar  to  Figures  8  through  11,  the  values  at  the  stopbands  are 
not  given.  For  better  shielding  effectiveness,  one  should  try  to  increase 
Z'Y,  and  avoid  md/c  =  nir  (n  =  1,2,3,  .  .  .  ). 


Figure  12.  Normalized  |z  T|  when  z  =d/2  versus  frequencies  for  various 


Figure  13.  Normalized  |«  |  when  z  =d/2  versus  frequencies  for  various 
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III.  DISTRIBUTED  EXCITATIONS 

In  this  section  both  the  outer  and  inner  shields  are  assumed  to  have 
uniformly  distributed  transfer  functions:  (shield  transfer  impedance 

per  unit  length) ,  (charge  transfer  frequency)  for  the  outer  shield 
and  Z^2>  for  the  inner  shield  (Figs.  2  and  14).  The  outer  shield  is 
different  from  that  discussed  in  Section  II  which  has  localized  transfer 
parameters.  The  total  current  I£  and  total  charge  per  unit  length  Q'  of 
the  double-shield  cable  are  assumed  known.  1^  and  Qj_  are  also  assumed  to 
be  dependent  upon  z  as  exp(-jhz)  and  related  to  each  other  via 

It  -  Jhlt  (27) 

In  what  follows,  results  for  this  distributed  double  shield  without  bonding 
straps  are  presented  first,  and  then,  the  more  general  case  with  periodic 
bonding  straps  is  discussed. 

1.  SCHELKUNOFF'S  CIRCUIT  (REF.  6) 

When  the  shields  are  solid  tubular  conductors,  =  and  the 

current-source  term  1^  of  the  transmission-line  equations  of  the  inside 

cables  vanishes.  As  for  the  voltage-source  term  V',  when  hd  =  0  and 

2  s 
o^iOi  >>  h  ,  it  can  be  calculated  from  the  circuit  depicted  in  Figure  15 

which  is  valid  regardless  whether  the  bondings  are  present  or  not 

(Ref.  6).  The  circuit  elements  in  the  figure  are  defined  and  given  as 

follows: 


Z*  =  surface  impedance  per  unit  length  of  the  i-th  shield 
ai 

with  internal  return  (i  =  l  for  the  outer  shield, 
i  ■  2  for  the  inner  shield) 

“  27.^-d;  [Io(nai)Kl(Yibi)  +  Ko(Yiai)Il(Yibi)] 
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Schelkunoff 's circuit  for  the  calculation  of  Z'  of  a  solid  tubular  double  shield 


*  surface  impedance  per  unit  length  of  the  i-th  shield 
with  external  return 


rbj-  [v-'iWvi’ +  Wi'WP 


“Ti  2jra.a.b.D, 
i  1  x  1 


=  series  impedance  per  unit  length  of  the  double 
shield  when  assumed  perfect 


=  jwpc|iln(a1/b2)/(2Tr) 


°1  ■ 


Y^  =  jwpiai  «■  2j/62 


pi’  °i’  3i»  bi»  =  Peri“eability,  conductivity,  inner 

radius,  outer  radius,  skin  depth 
of  the  i-th  shield 

and  I  ,  I, ,  K  ,  K  are  the  modified  Bessel  functions.  From  the  circuit, 
o  i  o  i 

one  immediately  finds  that  the  effective  transfer  impedance  per  unit  length 
ZT  Vs/V  the  shield  is  given  by 


where 


2n2;2 

*t  z;i+H2+z;i 


Z'  Z' 
T1  T2 


where 


zi  -  zk +  z;2  +  zk 


Evidently,  Equation  30  and  the  circuit  in  Figure  15  can  be  extended  easily 
to  describe  the  effective  transfer  impedance  per  unit  length  of  a  N-surface 
solid  tubular  shield  with  N  >  2. 
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Although  Equation  30  looks  simple,  the  circuit  elements  in  the  expression 
are  complicated  functions  of  frequency  and  shield  parameters  (see  Eqs.  28  and 
29).  When  A  (=  b.  -a.,  the  thickness  of  the  i-th  shield)  <<  a.,b.  and  Y.a., 

lii  ii  'll’ 

>>  1.  the  circuit  elements  Z^V,  Z^  ,  can  be  approximated  by 


z\  =  z'  =  z*  = 

ai  bi  Ii 


Vi 


2ir/a  b.a.A. 
i  i  l  i 


coth(Y.A.) 


i  l 


(32) 


z: 


Ti 


Y.A. 

- ^ -  csch(Y.A.) 

2*/Tb7o.A.  1  1 

1111 


(33) 


where  Z^  is  referred  to  as  the  internal  impedance  per  unit  length  of  the 
i-th  shield.  Equations  32  and  33  can  be  further  approximated  by 


Z*  =  Z'  -  Z'  -  Z'  =  R'  =  - - - 

ai  bi  ^Ii  ti  V,i  .  /— r-  . 

2ir/aibiaiAi 


(34) 


when  Y^A^  <<  1  (i.e.,  6^^  »  A^) ,  and 


Z' .  =  Z' ,=  Z’  =  Y.A.R'  j 
ai  bi  Ii  i  i  dc,i 


Zl.  «  2y.A.R'  .e 
Ti  'i  l  dc,i 


-y.a. 

1  i  a 


i  i 


(35) 


(36) 


when  y.A.  >>  1  (i.e.,  A.  >>  5.).  Here,  R'  .  is  the  dc-shield  resistance 
ii  ii  dc,i 

per  unit  length  of  the  i-th  shield  and  the  real  part  of  y^  is  taken  to  be 
positive. 

2.  CASEY'S  CIRCUIT  (REF.  7) 

When  <$.,  >>  A^  (which  guarantees  Z'  =  Z^.  =  Z|v  =  Z^ ,  see  Eq.  34), 

a  1 


"bi 


Ii 


Ti’ 


ry  *  A  A  A  WA  JL  -1.  1  A 

h  b  =  -  Z'Yjb,  -  -  ZVY'-b'  <<  1  (Y',  is  the  shunt  admittance  per  unit 
1  1  1  1  cl  cl  1  cl  r 

length  of  the  double  shield  when  assumed  perfect  and  no  bonding  straps),  and 
there  is  no  bonding  strap,  the  circuit  diagram  shown  in  Figure  16  can  be  used 
to  calculate  the  current  I  flowing  on  the  inner  shield  (Ref.  7).  In  the 
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circuit  (i  =  l,2)  are  the  capacitive  transfer  admittances  per  unit  length 
of  the  shields,  which  are  related  to  ft  via, 


where 


ft„.  =  -  sy:/y' 
Ti  i  Ti 


(37) 


(38) 


1 


Y' 

*2 


(39) 


and  Y^  is  the  shunt  admittance  per  unit  length  of  transmission  line  formed 
by  the  inner  shield  and  the  wires  within  when 

From  the  circuit,  one  has 


zh  *  h  /YT1 


<zti  +  h2/Yh>  +  <2n  +  h2/Yi2>  +  <zh  +  h2/Yu> 


where 


Z'  +  h  /Y' 
TI  TI 

zi  +  h2/Yi 


Zi  = 


z ' 
TI 


+  Z ' 
T2 


+  Z' 
cl 


from  which  one  immediately  has 


Z’  +  h' ft' 
TI _  Tl 

+  h2/Y^ 


(40) 


(41) 


(42) 


When  h  -*■  0  and  Y^  =  Y^  *  ®,  Equation  42,  indeed,  reduces  to  Equation  30 
with  6^  >>  A^.  It  is  believed  that  the  above  thin  shield  (6^  >>  A^) 
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approximation.  Equation  42,  can  also  be  used  for  thick  shields  provided  that 
Z|  defined  in  Equation  41  is  replaced  by  Equation  31,  i.e. , 


2i  =  Zk 


and  a  more  general  circuit  can  be  constructed  (Fig.  17).  Here,  the 

circuit  elements  Z'  and  Z'„  should  be  defined  in  a  broader  sense  than 
al  b2 

Equation  28  to  include  all  kinds  of  penetrations.  This  circuit  can  also 
be  easily  extended  to  describe  a  N-surface  shield  with  N  >  2. 

Up  to  this  point  the  discussion  in  Part  2  has  been  restricted  to  the 
situation  that  there  is  no  bonding  connecting  the  shields.  When  there  are 
periodic  bondings  with  period  d  «  h"1,  (/-YjzjT1,  it  is  postulated  that 
the  circuits  and  equations  can  still  be  used  provided  that  Y|  is  replaced 
by  Y^  +  Y^/d, i.e.,  Equation  42  becomes 


Z'  +  h  /Y* 

T  =  T1  Tl 

QT2  Z|  +  h2/(Y]^  +  Yd/d) 


(43) 


In  the  following,  a  general  analysis  of  a  double  shield  with  periodic  bondings 
will  be  given.  The  results  of  the  analysis  will  show  whether  the  simple 
Equation  43  is  accurate  enough  under  the  imposed  conditions. 

3.  GENERAL  FORMULATION 

To  obtain  the  current  I  flowing  on  the  inner  shield  one  has  to  solve 
the  following  transmission-line  equations  (Fig.  14): 


(44) 

(45) 
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The  shields  are  assumed  to  be  infinitely  extended  in  both  +z  and  -z 
directions  so  that  only  the  particular  solution  of  the  equations  need  be 
considered.  This  particular  solution  can  be  written  as  follows  (from 
Floquet ' s  theorem) 


Kz)  -  e‘jhz  l  ^e-j^/d 


V(z)  -  e'jhz  I  vPe“j2n7r2/d 


After  substituting  Equations  46  and  47  into  Equations  44  and  45,  and  going 
through  some  complicated  algebraic  manipulations  and  series  summations, 
one  eventually  has,  for  md  £  z  <  (m+l)d, 

Yi  1,2 * 2nYii  . ,  hYd  z«Yk-ziY; 

■*-Vz/  ^  y’  2  A  +  2Y1  ?X 

11  h  +2iYi  1  2nYii+h 

2  sin[(h+(-l)n/^zj")d/2]exp[-j  (h-(-l)n/=YTzT)  (md  -  z  +d/2)  ] 

x  l  -  1  1 - -  («) 

n-1  cos(d/-Yjzp  -  cos(hd)  +  YdZ|sin(d/-Yjzp / (2/-Y|zp 


from  which  Q'(z),  Z^,  and  can  be  calculated  via 


Q'(z) 


1  d  t  /  \ 
ju  dz  I(z) 


zi  =  z;2i(z)/it 


,  ,  J  T2  d_  I(z) 

T  h  dz  I 


Both  Z^  and  are  z-dependent.  In  order  to  better  quantify  the 
transfer  functions,  one  takes  the  average  values  of  Z^  and  over  the 
period  d  of  the  bonding  straps.  These  average  values,  designated  as  Z^ 
and  are  independent  of  z  and  are  given  as 


z;  $L  ipe^hz 
T  _  _T_  m  _ o _ 

ZT2  ^T2  Xt 


Y 1  h  +  Z '  Y ' 

1  ZT1  TI 

Y '  2 

TI  h  +  ZjY| 


2  Z'  Y'  -  7'Y’ 
1  .  h  Tl  T1  Z1Y1 

A  2  2 
n  +7'Yf  Z*  Y*  +  Vi 

11  zri  ti  ti 


cos  (d/-Y^Zp  -  cos(hd) 


dY,' 


1  cos(d/-Y[Z[)  -  cos(hd)  +  YdZjsin(d/-Yjzp/(2/^Yjz|) 


(52) 


It  can  be  shown  that  Equation  52  reduces  to  Equation  A3  when 


|  hd  -  /-Yp[  d 


lot  |  «  1,  hd  »  qs  («  Z^Ydd) 


(53) 


|  tan  hd  J  >>  a  ,  (i.e.,  hd  -  nirj  >>  |a[,  n  =  0,l,2,  ...  )  (5A) 


The  conditions  53  and  54,  obviously,  are  different  and  less  restrictive 

than  those  imposed  on  Equation  43  during  the  discussion  in  Part  2.  In 

most  practical  situations.  Conditions  53  and  54  can  be  satisfied,  and 

Equation  43  can  be  used.  In  the  case  that  the  constraints  53  and  54  are 

not  met,  one  has  to  resort  to  Equation  52  which  is  a  complicated  function 

of  the  frequency  and  the  s'vield  parameters.  The  shield  parameters  Z^, 

Y^i>  Y|,  Z|,  etc.,  are  generally  complex  values.  However,  when  the 

diffusion  penetration  is  not  important  (which  is  true  for  most  highly 

conducting  shields  at  frequencies  larger  than  10  kHz),  the  shield  parameters 

ZT1’  ^Tl’  etc,>  become  purely  imaginary  (i.e.,  the  penetration  is  through 

the  apertures  such  as  in  the  case  of  braided  cable  shields),  and  Equation  52 

becomes  a  real  function.  This  real  function  is  plotted  in  Figures  18  through 

25  as  a  function  of  hd  s  /-Y'  ,Z' .  d  *  wd/c  (from  0  to  10)  by  using  q  =Z', /Z',, 

cl  cl  z  n  ci 

qY  *  ^cl^Tl  Snd  *  Zcl^dd  0*5.  1,2)  as  parameters  (also  assume 

ZT2^Zcl  “  Ycl/,YT2  3  va^ues  f°r  the  parameters  (qY,qz)  *  (0.004, 

0.01)  and  (0.001,  0.002)  are  for  some  typical  braided  cable  shields  (see 
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100  =  0  frOOO 


Hgure  19.  Normalized  Z^,  and  (normalized  to  their  values  when  the  cable  shields  are 

not  bonded),  based  on  the  general  Equation  52  and  the  approximated  Equation  A3, 
versus  frequencies  for  qg  =  Z^Vjd-0.5,  qy  =  Y^/Y’ L  =  0.00A  and  q  =  Z '  /Z*  =0. 
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bonded),  based  on  the  general  Equation  52  and  the  approximated  Equation  43,  versus 
frequencies  for  q„ = Z'  Y  d = 2,  qv = Y'  /Y '  = 0. 004  and  q  = Z'  /Z'  = 0.01. 


Pages  581  and  584  of  Ref.  4).  In  the  figures,  the  curves  of  the  approximate 
Equation  43  are  also  given.  The  agreement  between  Equations  43  and  52  in 
the  region  where  wd/c  is  not  close  to  nn  (n=l,2,  .  .  .  )  is  clearly  shown 
even  when  hd  is  smaller  or  in  the  order  of  q^.  The  reason  they  agree  when 
hd  <  q  is  simply  because  both  Equations  43  and  52  have  values  in  the  order 
of  or  qz  (<<  1).  Because  of  this,  probably  one  can  lift  the  condition 

hd  >>  q  in  Equation  53.  From  the  figures  one  can  draw  the  same  conclusions 

J 

as  those  of  the  discrete  excitation  case  that  the  bondings  improve  the 

shielding  effectiveness  at  certain  frequency  ranges  while  degrade  it  at 

others.  In  order  to  widen  the  frequency  ranges  for  better  shielding 

effectiveness,  one  may  try  to  increase  Z'Y,.  Also  one  should  not  use  a 

1  d 

d  -value  which  causes  resonances  (where  tod/c  -  nir,  n  =  l,2,  .  .  .  )  and  thus 
seriously  degrades  the  shielding  effectiveness. 

Equation  52  can  also  be  rewritten  in  the  following  form: 


ZT1  *  "2/'1' 


T1 


ZT  “t 

7’  ~  q  ~  2 

i2  x2  Z[  +  h  / (Y^  +  qYd/d) 


where 


q  = 


1  +  dYdZ[(q2-qi)/[(h2  +  YjZpd2] 


qx  =  (Z^Y^  -  ZiYp/W^  +  h  ) 


(55) 


(56) 


(57) 


sin(d^-YJ Z' ) / (d/-Y'Z' ) 

q  =  - — ~ - 2^ - 2 -  (58) 

[cos(d/-Y^zp  -  cos(hd) ] / [ (n  +  YjZ|)d^/2] 


Equation  55  can  be  easily  represented  by  a  circuit  (Fig.  26).  Under  the 

assumption  that  the  diffusion  penetration  is  not  important,  the  q-values 

are  plotted  in  Figures  27  and  28  as  functions  of  hd  =  /-Y'  Z'  d  *  wd/c, 

cl  cl 

with  qz>  q^  and  qg  as  parameters.  The  values  of  the  parameters  in  Figures  27 


51 


double  shield  with 


100  0  = 


)02  and  various  q  =Z'  Y,d. 


and  28  are  the  same  as  those  of  Figures  18  through  25.  From  the  figures, 
one  can  easily  see  that  q  *1  when  Conditions  53  and  54  are  satisfied; 
that  is,  Equation  55  reduces  to  Equation  43.  This  is  another  proof  that 
the  simplified  Equation  43  can  be  used  for  Equation  52  or  55  when  the 
Conditions  53  and  54  are  satisfied. 
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IV.  SUMMARY 


In  this  section,  the  results  of  Sections  II  and  III  are  summarized 
and  some  examples  are  worked  out  to  show  how  the  results  can  be  used. 

1.  DISCRETE  EXCITATIONS 

When  the  outer  shield  of  a  double-shield  configuration  is  coupled 
to  a  localized  voltage  and/or  current  source  (Figs.  1  and  3),  one  may 
employ  the  Floquet  theorem  to  the  periodic  transmission-line  equations 
to  determine  the  disturbances  propagating  down  the  bonded  double  shield. 

A  passband-stopband  structure  in  the  dispersion  relation  between  w  and  k 
is  observed  (Eq.  8  and  Fig.  6).  In  the  stopbands  the  disturbances 
decay  exponentially  away  from  the  penetration  point,  whereas  in  the 
passbands  the  disturbances  oscillate  persistently.  In  the  stopband  the 
decaying  constant  can  be  easily  determined  from  the  dispersion  relation. 
Curves  of  the  decaying  constant  are  plotted  in  Figure  7.  At  low  frequencies 
the  decaying  constants  can  be  calculated  from  the  simple  approximate 
Equations  9  and  10. 

a.  Voltage  source  (Vq,  see  Figs.  3b  and  4a) 

The  voltage  source  can  be  calculated  from 


V 


s 

o 


IscZTl 


(59) 


where  I  is  the  short-circuit  current  on  the  outermost  surface  of  the 
sc 

double  shield,  and  ZT1  is  the  localized  transfer  impedance  of  the  outer 
shield.  This  voltage  source  gives  rise  to  V^  and  1^  which  are  the 
voltage-  and  current-sources  exciting  the  wires  inside  the  inner  shield. 
By  defining  a  combined  effective  transfer  impedance  per  unit  length  Z ^ 
and  a  combined  effective  charge  transfer  frequency  per  unit  length  as 


Z ' 
TV 


V'/I 
s  sc 


o' 

TV 


^uIs/,Isc 
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one  can  fully  describe  the  coupling  through  the  bonded  double  shield  due  to 

a  voltage  source  V®  (or  Igc>  see  Eq.  59).  and  are  averaged  over 

the  period  d  via  Equation  20.  The  averaged  quantities  Z^v  and  0^  are 

given  in  Equation  21  and  plotted  in  Figures  8  through  11  after  normalized 

-1/2 

to  their  values  for  no  bonding  strap  (ZtiZt2^Zi/YP  /np  anc*  ZTlfiT2YJ/np’ 
respectively,  where  np *  1  for  Figure  3b  and  np *  2  for  Figure  4,  and 
are  the  coupling  coefficients  of  the  inner  shield) . 

b.  Current  source  (Ia,  see  Figs.  3b  and  4b) 

The  current  source  can  be  calculated  from 


Q' 

sc  T1 


where  Qgc  is  the  short-circuit  charge  density  on  the  outermost  surface  of 

the  double  shield  and  0^  is  the  localized  charge  transfer  frequency  of  the 

outer  shield.  This  current  source  gives  rise  to  different  V'  and  I' 

s  s 

(different  from  those  due  to  V®).  By  defining  a  combined  effective  transfer 
impedance  Z^  and  a  combined  effective  charge  transfer  frequency  nTI  via 


ZTI  ' 


I'/Q' 
s  xsc 


one  can  describe  the  coupling  through  the  bonded  double  shield  due  to  a 

g 

current  source  Iq  (or  Qgc,  see  Eq.  60).  Z^  and  are  also  averaged  via 
Equation  20.  The  averaged  quantities  Z TI  and  are  given  in  Equation  26 
and  plotted  in  Figures  12  and  13  after  normalized  to  their  values  for  no 
bonding  strap  (flTjlZT2^np^“*  ’  nTiaT2^crV  resPectively>  • 

Figures  8  through  13  (or,  more  generally,  Eqs.  21  and  26)  show  that 
the  absolute  values  of  the  normalized  transfer  quantities  are  less  than  one 
(where  the  bonding  straps  improve  the  shielding  effectiveness)  at  some 
frequency  ranges  and  greater  than  one  at  the  others.  The  figures  also  show 
that  the  transfer  quantities  become  infinite  at  wd/c  =  nir  (n  -  1,2,3,  .  .  .). 
The  transfer  quantities,  however,  are  not  given  in  the  stopbands  (also  see 
Figs.  6  and  7)  where  they  are  extremely  small  far  away  from  the  penetration 


From  Che  above  observations,  one  concludes  that  in  order  to  have  a 
better  shielding  effectiveness,  the  periodic  bonding  should  be  employed 
such  that  the  important  part  of  the  EMP  spectrum  lies  inside  the  first 
stopband.  This  can  generally  be  realized  by  choosing  an  appropriate  d 
when  and  are  specified.  Also,  in  order  to  have  a  broader  first 
stopband  and  a  greater  decaying  constant,  one  should  try  to  make  ZjY^  larger. 

c.  Examples 

For  a  coaxial  double  shield  with  outer  radius  (b)  =  5  cm  and 
inner  radius  (a)  =  3  cm,  one  has 


,  h  V  7 

Z'l  =  ~ 2”  l)  "  10"  (0/») 


i.e. , 


L[  =  10“7H/m 


Also,  suppose  that  highly  conducting  wires  of  radius  t = 1  mm  are  to  be 
used  for  bonding,  then. 


Yd  ~  =  (ja))  M  (b-a)  Zn(2(b-a)/t) 

o 


i.e. , 


Thus, 


=  6.7  X  107(jw)"1  (r.)-1 


Ld  ®  1.5  x 10  h 


Z|Yd  =6.8  (m-1) 


The  question  now  arises  as  to  the  spacing  d  to  be  used.  Take  d  =  0.6m 
and  0.3m  as  given  in  Table  1.  Both  cases  give  rather  wide  first  stopbands 
which  cover  the  important  EMP  spectrum.  However,  the  case  of  d  *  0.3m 
gives  a  wider  first  stopband  and  a  larger  decaying  constant  (note  that  a 
decaying  constant  of  l(m  corresponds  to  an  atteunation  of  8.7  dB/m) . 
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Another  important  quantity  is  Z|Y^,  which  should  be  made  as  large  as 
possible  for  better  shielding  effectiveness.  If  the  double  shield  is 
coaxial,  one  can  use  a  smaller  (b  -  a)  to  obtain  a  larger  Z^Y^  (Eqs.  SI  and  62) 
The  case  of  b»4  cm  and  d  =  0.6  m  is  also  given  in  Table  1.  From  the 
table,  it  is  observed  that  the  first  stopband  is  wider  and  the  decaying 
constant  is  larger  for  b“4  cm,  d  =  0.6  m  than  the  case  for  b  *  5  cm,  d»0.6m. 


TABLE  1.  EXAMPLES  OF  BONDED  COAXIAL-CABLE  SHIELDS 


b  =  5  cm  a  = 

3  cm  t  =  1  mm 

b  =  4cm  a  =  3  cm  t  =  lmm 

(L*  =0.1yh/m,  L,  =  15  nh) 

1  a 

(L|  =  57  nh/m,  L^  =  6  nh) 

Period  of  Bondings 

d  =  0. 3  m 

d  =  0.6  m 

d  =  0.6  m 

%  <-  4/Ld> 

,  ,  -1 

6.7  m 

9.6  m-1 

1st  Stopband 

0-^205  MHz 

O-h-137  MHz 

O-m-160  MHz 

| Im(k) | ,  Below  10  MHz 

2.9  m-1 

3.4  m'1 

2.  DISTRIBUTED  EXCITATIONS 

Given  the  distributed  transfer  parameters  of  both  the  inner  and 
outer  shields  (Figs.  2  and  14),  the  effective  overall  transfer  parameters 
of  the  double  shield  can  be  represented  by  simple  circuit  diagrams.  These 
circuit  diagrams  are  summarized  below. 

a.  Schelkunoff ' s  circuit  (Ref.  6) 

When  the  shields  are  solid  tubular  conductors  whose  skin  depths, 
linear  cross-sectional  dimensions  and  the  period  of  the  bonding  straps 
are  much  smaller  than  the  wavelength  of  the  EMP  disturbance,  the  effective 
transfer  impedance  of  the  double  shield  (the  effective  charge  transfer 
frequency  0^,  is,  of  course,  zero)  is  independent  of  the  bonding  straps 
and  can  be  calculated  from  the  circuit  depicted  in  Figure  15  or  from 
Equation  30.  The  circuit  elements  in  Figure  15  and  Equation  30  are 
given  in  Equations  28  and  29,  and  also  Equations  32  through  36  for  some 
special  cases.  The  circuit  can  be  easily  extended  for  a  N-surface  solid 
tubular  shield  with  N  >  2. 


IP 


b.  Casey's  circuit  (Ref.  7) 

When  (a)  the  skin  depths  of  the  shields  are  much  greater  than 
their  thickness  (6  >>  A,  thin  shields),  (b)  the  linear  cross-sectional 
dimensions  of  the  double  shield  are  much  less  than  the  wavelength  of  the 
EMP  disturbance  (hb  =  /-Y^Z^  b  «  1) ,  and  (c)  there  is  no  bonding  strap, 
then,  Z^,  and  can  be  calculated  from  the  circuit  of  Figure  16  or  from 
Equation  42.  The  circuit  elements  in  the  figure  and  the  equation  are 
given  in  Equations  37,  38,  39, and  41. 

When  the  shields  are  not  thin,  the  circuit  of  Figure  17  (called 
the  generalized  Casey's  circuit)  can  be  used  to  replace  that  of  Figure  16 
for  calculating  Z^  and  Equation  42  is  still  applicable,  except  that 

the  circuit  element  Z^,  originally  given  by  Equation  41,  becomes 


7  * 

L 1 


z'  +•  z' 

al  d2 


z;i 


Here,  and  Z^2  are  defined  in  a  broader  sense  than  Equation  28  to  include 
all  kinds  of  penetrations.  Both  circuits  in  Figures  16  and  17  can  also 
be  easily  extended  to  describe  a  N-surface  shield  with  N  >  2. 

c.  General  circuit 

When  there  are  periodic  bondings  connecting  the  shields,  Z|  and 
(the  average  values  of  Z ^  and  QT  over  the  period  of  the  bondings)  can 
be  calculated  from  the  circuit  in  Figure  26  or  from  Equation  55  (or 
equivalently,  Equation  52).  The  constant  q  (qY^/d  may  be  named  the 
effective  shunt  admittance  per  unit  length  of  the  bondings)  in  Equation  55 
and  Figure  26  is  a  complicated  function  of  the  shield  and  bonding  para¬ 
meters  (Eqs.  56  through  58).  However,  when  the  conditions  given  in 
Equations  53  and  54  are  met,  q  is  approximately  equal  to  1,  i.e.. 

Equation  55  reduces  to  Equation  43.  The  truth  of  this  statement  is 
further  supported  by  the  curves  in  Figures  27  and  28  for  q  and  the 
curves  in  Figures  18  through  25  for  the  normalized  and  (normalized 
to  their  corresponding  values  for  no  bonding  strap,  i.e.,  Equation  55 
with  q  *  0) .  Actually,  in  Figures  18  through  25,  the  curves  of  Equation  43 
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agree  with  those  based  on  the  exact  Equation  55  even  in  the  region  where 
hd  <  q  (i.e.,  the  condition  hd  >>  q  in  Equation  53  is  violated,  and  q 
is  quite  different  from  one  as  shown  in  Figures  27  and  28),  provided  hd 
is  not  close  to  nit.  One  thus  concludes  that  when  q  >0.1,  ] hd  — 

/-YjZj  d  J<<1  and  hd  is  not  close  to  mr  (n  =  l,2,  .  .  .  )  Equation  43 
or  the  circuit  in  Figure  26  with  q  =  1  is  a  good  approximation  for 
calculating  Z^,  and  &T. 

From  the  results  presented,  one  concludes  that  in  order  to  shield 
against  the  distributed  excitations  more  effectively,  the  periodic 
bondings  should  be  implemented  in  such  a  way  that  Z'^Y^  is  large  and 
wd/c  <  tt  for  the  important  parts  of  the  EMP  spectrum. 

d.  Examples 

Consider  a  coaxial  double  shield  (Fig.  2)  with  outer  radius 

b  =  5  cm,  inner  radius  a  =  3  cm  and  with  highly  conducting  bonding  straps 

of  radius  t  =  1  mm.  Also,  take  the  practical  values  q^  =  Y^/Y^  =  0.004 

and  q„  =  Z'/Z*  -  0.01.  Then, 

Z  T1  cl 

Z’1Yd  -  6.7  (nf1) 

from  which 

2  for  d  =  0. 3  m 
1  for  d  =  0.15  m 

From  Figures  20  and  21,  one  immediately  sees  that  the  normalized  transfer 
functions  are  less  than  0.1  (i.e.,  the  bonding  straps  reduce  the  EMP 
penetration  by  more  than  20  dB)  for 


Z' -Yd  = 
cl  d 


frequencies  < 


470  MHz 
940  MHz 


when  d  =  0 . 3  m 
when  d  *  0. 15  m 


Both  470  MHz  and  940  MHz  are  extremely  high  to  be  important. 
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However,  tf  the  bonding  straps  between  the  coaxial  double  shield  give 
a  much  smaller  value,  say, 

Z ' , Y  .  =  0.67  m-1 
cl  d 

then,  when  d  =  0.15m  is  used  the  same  quality  in  the  shielding  effective¬ 
ness  as  that  of  Z^Y^  =  6.7  m-'*'  and  d  =  0.15  m  can  be  obtained  only  for 
frequencies  up  to  at  most  300  MHz  (see  Fig.  18).  Thus,  bonding  straps 
with  greater  Z^  Yd  are  preferred. 
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I.  INTRODUCTION 


One  of  the  most  effective  methods  of  protecting  ^  system  from  undesirable 
EMP  effects  is  by  use  of  shielding.  Although  the  basic  principles  of 
shielding  for  a  single-surface  enclosure  are  well  understood  (Ref.  1),  the 
extension  of  these  principles  to  multisurface  enclosures  is  by  no  means 
straightforward.  Unique  to  multisurface  shielding  is  the  mutual  interaction 
among  the  shields,  which  may  degrade  the  intended  shielding  performance  of 
the  enclosure.  Another  unique  feature  is  the  bonding  that  is  often  employed 
between  the  shields  in  order  to  reduce  electrostatic  hazards.  This  bonding 
practice  may  have  an  adverse  effect  on  the  protection  of  a  shielded  enclosure 
against  magnetic-field  penetration.  The  effects  of  shield-shield  interaction 
and  bonding  will  be  treated  in  this  report. 

Figure  1  shows  various  topics  that  will  be  addressed  in  this  report. 

In  Section  II,  the  problem  of  two  concentric  spherical  shields  will  be  solved 
using  the  theory  of  inductive  shielding  (Ref.l),  and  the  results  will  be 
generalized  to  N-surface  spherical  shields.  The  corresponding  results  for 
cylindrical  shields  will  be  presented  in  Section  III.  Equivalent  circuits 
will  be  constructed  in  Section  IV  to  interpret  the  analytical  results  for 
a  two-surface  spherical  enclosure;  the  results  will  be  generalized  to  two- 
surface  enclosures  of  arbitrary  shape.  In  Section  VI  the  effect  of  bonding 
on  magnetic-field  penetration  into  a  two-surface  enclosure  will  be  discussed. 
Finally,  the  most  important  results  are  summarized  in  Section  VII. 

The  underlying  assumptions  of  inductive  shielding  are  (1)  the  electric 
field  is  neglected  everywhere  except  in  the  enclosure's  wall  where  it  is 
related  to  the  induced  current  by  Ohm's  law,  and  (2)  the  wall  thickness 
is  much  smaller  than  the  typical  linear  dimension  of  the  enclosure.  In 
addition,  this  report  assumes  the  wall  thickness  to  be  smaller  than  the 
wall's  skin  depth,  except  in  Section  V  where  this  assumption  is  removed. 
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II .  SPHERICAL  SHIELDS 


In  this  section,  the  problem  of  two  concentric  spherical  shells  will 
first  be  treated  with  explicit  simple  engineering  results  given  in  frequency 
and  time  domains.  The  results  for  two  shells  will  then  be  generalized  to  N 
concentric  spherical  shells.  Discussions  of  equivalent  circuits,  generalization 
to  shields  of  arbitrary  shape,  and  effects  of  electrical  bonding  between  the 
shells  will  be  relegated  to  later  sections. 

1.  TWO  SPHERICAL  SHIELDS 

Figures  2a  and  2b  show  an  enclosure  with  two  concentric  spherical  shields 
immersed  in  a  slowly  varying  magnetic  field  HQ(t).  Insofar  as  the  penetrant 
field  H^(t)  is  concerned,  one  may  replace  Figure  2a  with  Figure  2b  with 
appropriate  boundary  conditions  that  duplicate  the  shielding  properties  of 
the  walls  (Ref.  1).  The  magnetic  scalar  potential  $  for  the  three  regions  shown 
in  Figure  2b  takes  the  form 

3 

'  al 

=>  -  H  r  cos0  +  A  —r  cos0  r  >  a, 

-L  O  ^4  1 

3 

ai 

■»  B  r  cos0  +  C  —  cos0  r  i  a2 

r 

*2  “  D  r  cos0  r  a^  (1) 

where  Hq  is  the  Laplace  transform  of  HQ(t).  The  constants  A,  B,  C  and  D  are 
determined  by  the  following  boundary  conditions  (Ref.  1  and  Appendix  A): 


3_ 

3r 


3r  *2’ 


at  r  -  a^ 


2 


i_  * 

3r  3’ 


at  r  -  a2 


(2) 
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7S  <W  "  SVl°l  h  *1*  at  r-aL 

7s  (W  “  sV2a2  f?  *2*  at  r-a2  (3) 

Equations  2  mean  that  the  normal  component  of  the  magnetic  field  is  continuous 
across  the  shield,  while  Equations  3  state  that  its  tangential  component  is 
discontinuous  by  the  amount  of  the  current  induced  in  the  shield.  From 
Equations  1  through  3  one  finds  that 


with 


HQ  (1+t1s)(1  +  t2s)  -  (a2/a^)  ^t1t2s^ 

T1  =  1  VlVl 

T2  *  1  voa2°2^2 


(4) 


(5) 


If  no  interaction  between  the  shields  is  assumed.  Equation  becomes 


(1  +  t,s)(1  +  t~s) 

no  int 

as  one  would  expect,  since  Equation  6  is  the  product  of  the  transfer  function 
of  each  individual  shield.  Equations  4  and  6  are  plotted  in  Figure  3. 

The  time-domain  solution  of  Equation  4  is 

*i<*> _ 1 

H°  /(  ^  2  ,  3 

y  (xx  -  x2)  +  ^TiT2a 
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Frequency  spectrum  of  the  ratio 


where 


2^(1- a  ) 


1»2 


rl+  t2  +  *^Ti“  T; 


2  3 

)  +  4T^T2a 


and  a  =  a2/a^,  while  the  time-domain  solution  of  Equation  6  is 


Ht(t) 

~H 


no  int 


(8) 


(9) 


Here,  Hq  is  the  impulse  strength  of  the  external  fields,  and  for  most 
applications  can  be  taken  to  be  the  time-integral  of  the  magnetic  field  of 
a  typical  high-altitude  EMP  (Ref.  1). 

Equations  7  and  9  are  plotted  in  Figure  4  where  one  may  see  that  for  the 
case  a2/a^  ~  neglect  of  shield-shield  interaction  amounts  to  20% 

underestimate  of  the  penetrant  field. 

From  the  viewpoint  of  the  EMP  hardness  designer  the  currents  induced 
in  each  enclosure's  shield  are  important,  since  they  are  the  only  means  to 
prevent  the  external  field  from  penetrating  into  the  interior  of  the  enclosure. 
Let  and  K2<,  denote  the  induced  sheet  currents  in  first  and  second  shields 
of  Figure  2.  Then,  from  Equations  1  through  3  one  obtains 


7  .  1  3$2  1  3*1 

Kl<p  ~  a1  30  “  a1  38  (at  r  ~  al) 


■  -  2®o 


3  2 

T1s(l+T2s)-a  T^T2S 
(1 + T^s) (1 + T2s)  -  a^T^T2S^ 


sine 


(10) 


K  ^  1  3*3  1  H2  (  . 

K2$  a2  36  “  a2  39  (at  r * a2) 


-  -  |  H 
2  o 


t2S 


3  2 

(1 + t^s) (1 + t2s)  -  a  t^s 


sin9 


(ID 
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field  for  an  impulse 


These  expressions  for  the  currents  will  give  a  clue  to  construct  equivalent 
circuits  in  Section  IV. 

2.  N  SPHERICAL  SHIELDS 

Let  H^ ^ ,  H^3^ ,  ...  H*N)  denote,  respectively,  the  field  that  penetrates 
into  a  two-surface,  three-SUrface >  _  N-surface  spherical  enclosure. 

For  a  two-surface  enclosure  Equation  4  gives 

so/s<2>  -  <l+V)U+t2S)  TlV2 

*  1  +  (t1  +  t2)s  +  [l-  (a2/a1)3  ]t1t2s2  (12) 

For  a  three- surface  shielded  enclosure  (Fig.  5a)  one  uses  the  same, 
although  more  complicated,  procedure  for  the  two- surface  enclosure  and  finds 

**o^i3)  *  1  +  ^Tl+T2  +  T3^s  +[1  “  (a2/ai>3  ]TlT2s2+  [1-  Ca3/a2)  3Jt2T3s2 

+  [  1  ~  (a3/ai)  3Jt3tis2  +  ~  Ca2^al^  3  J  [-1-  "  ^a3/,a2^3]TlT2T3s2 

The  poles  of  H^(s)  in  the  complex  s-plane  for  two  and  three  spherical 

shields  are  given  in  Table  1.  These  poles  will  immediately  enable  one  to 

plot  the  frequency  spectrum  of  the  penetrant  field,  since  they  are  the 

"break  points"  in  the  log-log  scale  plot  (Fig.  6).  From  the  table  it  is 

clear  that  as  the  second  and/or  the  third  shield  get  closer  to  the  outer 

shield,  the  pole  corresponding  to  the  outermost  shield  moves  toward  the  jw- 

axis  away  from  its  unperturbed  value  -1,  while  the  pole  corresponding  to  the 

second  (third)  shield  (the  second  (third)  column  of  Table  1)  moves  away  from 

2 

its  unperturbed  value  -1/a  (-1/a  )  further  away  from  the  jw-axis. 

For  an  N-surface  spherical  enclosure  (Fig.  5b)  one  can  write  down,  on 
a  close  examination  of  Equations  12  and  13, 
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j(cu)/H0(cu)l 


TABLE  1.  s- PLANE  POLES  FOR  TWO- SURFACE  AND  THREE- SURFACE  SPHERICAL 
ENCLOSURES.  THE  VALUES  ARE  FOR  s  =  is  WITH 
a  =  a„/a„  =  a,/a„  AND  x„/x.  =  x,/x„  =  a 


Two 

shields 

Three  Shields 

a 

?1 

¥2 

S1 

s2 

S3 

0.1 

-  1.00 

-10.01 

-  1.00 

-10.01 

-100.11 

0.2 

-  1.00 

-  5.05 

-  1.00 

-  5.04 

-  25.25 

0.3 

-  0.99 

-  3.46 

-  0.99 

-  3.43 

-  11.55 

0.4 

-  0.96 

-  2.77 

-  0.96 

-  2.67 

-  6.96 

0.5 

-  0.91 

-  2.52 

-  0.90 

-  2.29 

-  5.10 

0.6 

-  0.83 

-  2.57 

-  0.79 

-  2.13 

-  4.49 

0.7 

-  0.73 

-  2.96 

-  0.66 

-  2.17 

-  4.72 

0.8 

-  0.65 

-  3.96 

-  0.53 

-  2.56 

-  6.03 

0.9 

-  0.57 

-  7.22 

-  0.42 

-  4.10 

-  10.81 

aj  =  air 


Figure  6.  Frequency  spectrum  asymptotes  and  break  points  where  to 
&)-,  w-  are  given  in  Table  1  (ST,  =  Is,/,  SL  =  IF,  I,  etc 


N  2  N,N 


1  +  s  £  t  +  s 

i-1  i> j 


1  -  (a^j^/a  )' 


T  .  T  .  + 

1  J 


N 


n  r 


_ +  s"  n 

i»l 


1  “  (ai/ai-ir 


T  . 
1 


(14) 


*  (N) 

where  aQ  =  As  expected,  has  N  poles  lying  on  the  negative  real  axis 

of  the  s-plane. 
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III.  CYLINDRICAL  SHIELDS 


The  procedure  of  solving  the  problem  of  a  multisurface  cylindrical  shield 
where  the  external  magnetic  field  is  perpendicular  to  the  axis  of  the  shield 
follows  exactly  that  of  the  spherical  shield  described  in  Section  II.  The 
other  polarization  where  the  external  magnetic  field  is  parallel  to  the  axis  of 
the  shield  is  treated  in  Reference  2.  From  the  geometry  depicted  in  Figures 
7a, d  one  can  immediately  write  down 


$ 

$ 


1 

2 


H  p  cosh  +  —  cosh 

o  p 


=  Bp  cos4>  ■+•  —  cos<p 


P  L  bl 
bl  —  P  —  b2 


(15) 


$2  =  Dp  cos<j> 


P  1  b2 


Applying  at  p  *  b^  and  p  =  b2  the  boundary  conditions  (Ref.l) 


3  .  3  . 

■  -  <T>  »  ■  ■  © 

3p  1  3p  2 


P  =  b, 


—  $  =  —  $ 
3p  2  3p  3 


P  “  br 


^<*2  -  ■  sti  k  *1 


P  =  b. 


Vs<‘3  -  V  '  St2  k  ‘ 2 


p  =  b. 


(16) 


to  Equations  15  one  gets 


Hi(s)  L 

—  ’ ■  —  »—  3  ■  ■  -■  ■  ■  ■  —  ■  ~  ~  ~  I 

Ha(s)  (l+st1)  (l+st2)  -  (b2/b1)  2t1t2sii 


(17) 


where 
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Figure  7.  Muitisurface  cylindrical  shields. 


C1  =  2  Uobl°lAl 


C2  =  2  Mob2a2A2 


Similarly,  for  a  three-surface  cylindrical  shield  (Fig.  7b)  one  has 

Ho(s)/H^3)(s)  -  1+  (tL+t2+t3)s  +  [l-  (b2/b1)2]t1t2s2+  [l-  (b3/b2)2]t2t3s' 
+  [l  -  (b3/b1)2  j t3tLs2  +  [l  -  (b2/b1)2][l  -  (b3/b2)2]t1t2t; 

a 

and  for  an  N-surface  cylindrical  shield  (Fig.  7c)  one  has 


+  .... 


Hq(s)/H<N)(s)  =  1  +  s  ^  ti  +  s2  Y  [l  -  (b./b.)2]  t.t. 

....+s8ini[1-(b1/bi_1)2]ti 


where  bQ  *  A  comparison  of  Equations  17  through  20  with  Equations  4,  5,  13, 
and  14  reveals  that  the  results  for  spherical  shields  obtained  in  the  last 
section  can  be  directly  used  for  cylindrical  shields  if  one  replaces 


spherical  shields 


(aj/a^)' 


cylindrical  shields 


<VV 


Therefore,  in  the  following  two  sections  on  equivalent-circuit  representation 
and  effects  of  bonding, discussions  will  be  restricted  only  to  the  case  of 
spherical  shields. 

Table  2  gives  the  s-plane  poles  for  for  two  and  three  cylindrical 
shields,  while  Figure  8  shows  the  frequency  spectrum  asymptotes  and  break 
points  for  |H^/H  | . 
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TABLE  2.  s-PLANE  POLES  FOR  TWO-SURFACE  AND  THREE- SURFACE  CYLINDRICAL 


SHIELDS.  THE  VALUES  ARE  FOR  s  =  st1  WITH  B  =  ^2^1* 
b.j/b2  AND  t2/t  =  t^/t2  *  B 


B 

Two  Shields 

Three  Shields 

S1 

s2 

S1 

s2 

s3 

0.1 

-  1.00 

-10.11 

-  1.00 

-10.10 

-101.12 

0.2 

-  0.99 

-  5.26 

-  0.99 

-  5.21 

-  26.30 

0.3 

-  0.96 

-  3.80 

-  0.96 

-  3.66 

-  12.68 

0.4 

-  0.92 

-  3.25 

-  0.91 

-  2.98 

-  8.20 

0.5 

-  0.85 

-  3.15 

-  0.82 

-  !.67 

-  6.51 

0.6 

-  0.77 

-  3.40 

-  0.71 

-  2.60 

-  6.13 

0.7 

-  0.69 

-  4.07 

-  0.59 

-  2.30 

-  6.74 

0.8 

-  0.62 

-  5.63 

-  0.49 

-  3.47 

-  8.85 

0.9 

-  0.55 

-10.56 

-  0.40 

-  5.85 

-  16.09 

.t. 


IV.  EQUIVALENT  CIRCUIT  REPRESENTATION  -  GENERALIZATION 
TO  SHIELDS  OF  ARBITRARY  SHAPE 

The  results  obtained  in  the  last  two  sections  will  be  interpreted  in 
terms  of  equivalent  circuits  in  this  section.  The  advantages  of  equivalent 
circuit  representation  of  mathematical  results  are  two-fold:  (1)  it  is  useful 
for  interpreting  results  and  understanding  physical  mechanism  involved,  and 
(2)  it  is  a  quick  way  to  generalize  the  results  for  specific  shapes  of  en¬ 
closure  to  arbitrary  shapes  of  enclosure. 

To  gain  more  familiarity  with  what  follows  one  starts  with  one-surface 
spherical  shielded  enclosure  (Fig.  9a)  whose  low-frequency  transfer  function 
is  (Ref .  1) 


with 


H^s) 

Ho(a) 


1 

1+TS 


k  W  acA 

3  o 


(21) 


(22) 


where  a,  a  and  A  are  defined  in  Figure  9a.  Equation  21  can  be  represented 
by  either  the  equivalent  circuit  of  Figure  9b  or  Figure  9c  where 


L 


3  V’ 


(23) 


The  inductance  L  can  be  expressed  in  terms  of  the  volume  V  and  the  surface  S 
of  the  enclosure  as 


L  ■  p  V/S  (24) 

o 

which  also  applies  directly  to  cylindrical  as  well  as  two-parallel-plate 
enclosures  (Ref.l).  Equation  21  can  thus  be  used  as  the  transfer  function 
of  a  single-surface  shielded  enclosure  of  arbitrary  shape  if  t  is  interpreted 
as 
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(25) 


with  L  given  by  Equation  24  and  R  by  (aA) 

The  extension  of  the  equivalent  circuit  given  in  either  Figure  9b  or  9c 
for  a  single-surface  enclosure  to  a  two- surface  enclosure  turns  out  to  be  a 
nontrivial  matter.  In  principle,  one  may  start  with  the  transfer  function 
given  in  Equation  4  and  constructs  an  equivalent  circuit  for  it  using  the 
techniques  known  in  circuit  synthesis  (Ref.  3).  This  approach,  however,  does 
not  easily  lead  to  a  circuit  which  represents  the  actual  physical  phenomenon 
of  the  problem.  To  derive  the  desirable  circuit  one  returns  to  the  induced 
currents  in  the  shields  given  by  Equations  10  and  11.  The  total  induced 
current  in  each  shield  is  obtained  by  integrating  Equations  10  and  11.  Thus, 


K,  a.de  =  -  3H  a, 
lq>  1  o  1 


3  2 

t^s(1 + x2s)  -  a  t^t2s 
(1  +  t1s)(1  +  t2s)  -  (a2/a1)3T1x2s2 


(26) 


-ft 


K„,a„d9  =  -  3H  a 


V 


j  2<t  2 


o  2  .  W1  ^  .  3  2 

(1 + t^s) (1 + t2s)  -  a  t^t2s 


(27) 


To  generalize  Equations  26  and  27  to  two- surface  enclosures  of  arbitrary 
shape  one  simply  sets 


T1 


a3T1T2  =  M2/(R1R2) 


(28) 


where  the  self-inductances  L^,L2,  the  mutual  inductance  M,and  the  resistances 
R^,R2  are  given  by 


Vi 


M 


Vx  L1L2 


Vl 


°2A2 


(29) 


Here,  and  are  the  volumes  enclosed  respectively  by  the  surfaces  and  S2 
of  the  first  and  second  shield.  Equations  26  and  27  can  be  written  in  the 
general  form 

sL.  (sL,  +  RJ  -  s2M2 

l - i - £ - £ -  !  (30) 

(sL1  +  R1)(sL2  +  R2)  -  s  M 

sMR- 

i2  - - - - r-r  i  (3i) 

(sL1+ R1)(sL2  +  R2)  -  s  M 

with 

Io  -  -  38^  (32) 

It  can  be  easily  verified  that  the  equivalent  circuit  shown  in  Figure  10  leads 
to  Equations  30  and  31. 

It  remains  to  show  how  1^  and  I2  are  related  to  the  penetrant  field  iL 
given  by  Equation  4,  which  can  be  expressed  in  the  generalized  form 


R1R2 


2..2 


Hq  (sL1  +  R;i)(sL2  +  R0)  -  s  M' 


(33) 


From  Equations  30  and  31  one  gets 


The  scaling  factor  L2/M  can  be  expressed  in  terms  of  the  geometric  parameters 
of  the  two  shields  with  Equations  29  and  is  given  by 


\  _  _ R1R2 _ 

/  (sL1  +  R1)(sL2  +  R2)  -  s2M2 


(35) 
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Hence,  one  has  on  comparing  Equations  33  and  34 


(36) 


which  is  to  say,  one  can  first  calculate  1^,  and  from  the  equivalent  circuit 
of  Figure  10  from  which  the  penetrant  field  is  directly  deduced  from 
Equation  36. 

A  final  point  should  be  made  about  the  scaling  factor  /sTJTs^  in  Equation 
36.  This  factor  comes  about  because  the  magnetic  field  H  is  proportional  to 
current  density  K  rather  than  the  total  current  I. 
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V.  TWO  SPHERICAL  SHIELDS  OF  ARBITRARY  ELECTRICAL  THICKNESS 


In  the  previous  sections  the  thickness  of  the  shield's  wall  was  assumed 
to  be  electrically  thin.  This  assumption  holds  true  for  low  frequencies 
and/or  poorly  conducting  shields.  In  this  section  this  assumption  will  be 
removed  and  the  shield's  wall  can  be  of  arbitrary  electrical  thickness. 

As  in  Equation  1  let  the  scalar  potential  $  for  the  three  regions  shown 
in  Figure  2b  take  the  form 

3 

al 

=  -  H  r  cos0  +  A'  — =•  cos0  r  >  a, 

1  o  2  —  1 

r 

*1 

$2  =  B'r  cos0  +  C'  — j  cos0  r  L  a2 

r 


$2  =  D'r  cos0 


r  <  a„ 


(37) 


Instead  of  the  boundary  conditions  given  by  Equations  2  and  3,  the  boundary 
conditions  are  now  given  by  (Ref.  1) 


3_ 

3r 

3_ 

3r 

3_ 

3r 

3_ 

3r 

where 


($2  +  ^)  =  aj7‘(*2  -  *j) 
($2  -  =  31V^($2  +  $1) 

($2  +  <f2)  =  a27s^$3  "  $2^ 

(*3  “  *2)  =  62Vg($3  +  $2) 

»  1 
ai  woPi  tanh(Pj/2) 

>JiAi 

“  — —  tanh(p  /2) 

1  uopi  1 

2  . 2 
Pi  -  suiaiAi 


at  r  =  a-^ 

at  r  =  (38) 

at  r  * a2 

at  r  -  a2  (39) 

(40) 
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and  i-1,2.  The  constants  A' , 
Equation  37  into  Equations  38, 
penetrant  field  is  obtained 


B',  C'  and  D'  can  be  found  by  substituting 
39  and  40,  and  are  given  in  Appendix  A.  The 
from  and  given  by 


H0  (cosh  Pj^  +  K1p1sinh  p^  (cosh  p2  +K2p2sinh  p2>-  (a^a^  KjPjK^sinh  p^sinh  p2 

(41) 


with 


K,  = 


“oal 


K„  = 


Uoa2 


1  ’  2  3lJ2A2 


(42) 


Just  like  Equation  4  the  combined  transfer  function  of  two  shields  is  the 
reciprocal  of  the  product  of  the  transfer  functions  of  individual  shields 
minus  an  interaction  term. 


For  electrically  thin  shield  walls  (i.e.,  p^,p2  <<:  l)*  Equation  41 
reduces  to  Equation  4,  as  it  should.  On  the  other  hand,  if  the  shield's 
walls  are  electrically  thick  (i.e.,  P1»P2  >>  1) »  then  Equation  41  gives 


H . 

l 


4e 


-(P1+P2) 


(l+KlPl)(l+K2P2)  -  (a2/a1)JK1P1K2p2 


(43) 


in  which  one  may  drop  the  ones  in  comparison  with  K^p^  and  K2p2>  However, 
Equation  43  is  preferable  since  it  is  similar  in  form  to  Equation  4  for  the 
case  of  an  electrically  thin  shield. 

Equation  43  is  plotted  in  Figure  11  with  and  without  the  interaction 
3 

term  (a2/a^)  K^p^K2p2>  It  can  be  concluded  that  for  a2/a^  =  0.9  and  ojt2 > 100 
(which  corresponds  to  f  >  16  kHz  for  t2  -  1  ms) ,  the  shield-shield  inter¬ 
action  reduces  the  shielding  effectiveness  by  at  least  a  factor  of  3  or  10  dB. 
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VI.  EFFECT  OF  BONDING 


AlChough  bonding  is  often  employed  between  shields  in  engineering 
practice  to  reduce  electrostatic  hazards,  there  is  no  quantitative  information 
on  its  possible  effect  on  the  shielding  performance  of  a  shielded  enclosure 
against  the  low-frequency  magnetic-field  penetration.  This  section  is 
devoted  to  calculating  this  effect  for  two  electrically-thin  spherical 
shields . 

Figure  12  shows  various  types  of  connection  arrangements  of  bonding 
straps  between  the  two  surfaces  of  a  spherical  enclosure.  Later  in  this 
section  it  will  become  clear  that  the  bonding  strap  arrangements  (a)  through 
(d)  have  no  effect  on  the  magnetic-field  shielding  performance  of  the 
enclosure,  while  the  arrangement  (e)  has  a  significant  adverse  effect. 

Before  proceeding  it  is  appropriate  to  remark,  that  the  low-frequency  electric 
field  within  a  conducting  enclosure  comes  mainly  from  the  eddy  currents 
in  the  enclosure's  wall  induced  by  the  time  rate  of  change  of  the  external 
magnetic  field.  It  is  this  electric  field  that  is  affected  by  the  bonding 
straps  in  the  inductive  shielding  approximation. 

Return  now  to  Figure  12e  (which  is  redrawn  in  Figure  13)  and  calculate 
the  current  induced  in  the  bonding  straps.  Let  the  unprimed  quantities 
be  the  quantities  in  the  absence  of  the  bonding  straps,  and  the  primed 
quantities  the  quantities  due  to  the  presence  of  the  bonding  straps.  Then, 
integrating  the  equation 


V  x  E  =  -  suq5  (44) 

over  the  area  enclosed  and  traced  out  counter  clockwise  by  the  loop  BCDAB 
of  Figure  13,  one  gets 


u  +  v'  -  V  -  V' 
BC  BC  AD  AD 


SUI  +  sv< 

ABCDA  ABCDA 


(45) 


where  V  is  the  voltage  drop  and  is  the  magnetic  flux.  In  deriving 
Equation  45  the  bonding  straps  have  been  assumed  to  be  good  conductors; 
otherwise,  a  term  for  the  voltage  drop  along  AB  and  CD  has  to  be  added  to 
the  left  hand  side  of  the  equation. 
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bonding  strop 


A  two— spherical  shield  enclosure  with  two  bonding  sti 
plane  perpendicular  to  the  external  magnetic  field  H 


Let  1^  ■  current  in  bonding  strap,  ■  inductance  of  one  bonding  strap, 
*BC  or 

shield.  It  is  obvious  from  Equation  44  that 


Rda  *  resistance  between  B,C  of  the  outer  shield,  or  D,A  of  the  inner 


^b(RBC+RDA  +  2sLb)  *  Sli,ABCDA  "  VBC  +  VAD 


oVU  _  qW  4-  ell/ 

s  ABCDA  OBCO  s  OADO 


“  s(¥OAB  +  ^OCD^ 


-  2s4< 


OAB 


(46) 


The  second  step  of  the  right  hand  side  follows  from  repeated  applications  of 
Equation  44.  It  has  been  assumed  for  simplicity  that  the  two  bonding  straps 
are  identical.  Hence,  the  triangle  OAB  is  equivalent  to  the  triangle  OCD. 

Similarly,  an  application  of  Equation  44  to  the  area  enclosed  by  the 
loop  BADCB  going  counter  clockwise  in  Figure  13  gives  Equation  46,  as  it 
should . 

Solving  Equation  46  for  1^  and  using  the  expression  in  Appendix  B  for 

„  one  obtains 
OAd 

i _ -.i!V 

b  rbc  +  rm+2sLi, 


FH  s(l  +  sT  ) 

o  o _ 

Rg  (1  +  sTfa) (1 +  sT1> (1  +  sT2) 


(47) 


where  ^b“2Lb^Rs’  T1  anc*  ^2  are  8iven  in  Squat*-00  8,  F  and  Tq 

are  given  in  Appendix  B.  It  can  be  shown  that  for  practical  cases  Tq  is 
small  compared  to  the  decay  time  constant  t2  of  the  inner  shield.  Hence 
for  low  frequencies  one  may  use  the  approximation 


ST. 


(l  +  sTb)(l  +  sT1)(l  +  sT2) 


(48) 
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with  VQ  »  FHq/t2»  where  FHQ  is  identically  equal  to  the  external  magnetic 
flux  linking  the  two  triangles  shown  in  Figure  13. 

The  magnitude  of  Equation  48  is  plotted  in  Figure  14  against  with 

Tt/T2  as  a  parameter.  The  time  histories  of  1^  and  1^  are  shown  in 
Figures  15  and  16  for  an  impulsive  external  magnetic  field  HQS(t),  which 
is  a  valid  representation  of  any  pulsed  external  field  whose  pulse  width  is 
less  than  the  diffusion  time  through  the  shield's  wall.  The  diffusion  time 
of  a  typical  metallic  enclosure  is  on  the  order  of  tens  of  microseconds.  The 
parameter  Tb/t2  is  roughly  equal  to  the  ratio  of  the  inductance  of  the  two 
bonding  straps  to  the  inductance  of  the  inner  enclosure,  the  latter  being 
given  by  pQa 2/3;  that  is  to  say 


T  2L 
b  >  b 


6Lb 

uoa2 


(49) 


where  can  be  estimated  from  the  approximate  formula 


Lb 


2rr 


£n(i/r) 


(50) 


with  £  *  length  of  one  bonding  strap,  r  ■  effective  cross-sectional  radius 
of  the  strap.  It  can  be  seen  from  Equations  49  and  50  that  Tb/r2  is  usually 
less  than  unity.  The  smaller  is  this  parameter  the  more  current  will  be 
induced  in  the  bonding  straps,  as  can  be  observed  in  Figures  14  through  16. 
Of  course,  the  more  current  there  is  in  the  bonding  straps  the  more 
penetration  there  is  into  the  enclosure.  Table  3  summarizes  the  peak 
values  for  Ife(t)  and  Ife(t)  for  various  values  of  Tb/r2.  In  the  table  it 
is  assumed  that  a2/a^«0.9,  Rg  a  (oA)  \  and  the  two  shields  have  the  same 
thickness,  conductivity  and  permeability. 

To  get  some  rough  estimate  on  the  field  due  to  the  bonding  strap 
current  Ife  one  may  divide  Ib  by  2ira2.  From  Table  3  and  for  the  case 
Tb^T2  “  t*ie  Peak  penetrant  fields  due  to  Ifa  are 


Figure  15.  Time  behavior  of  the  bonding  strap  current  due  to  an  impulsive  external 

2 

magnetic  field  Ho5(t).  The  normalizing  factor  is  VQ/Rg  with  V^FlW^ 
Note  that  H  has  the  dimension  of  ampere-second. 


Time  behavior  of  the  time  rate  of  change  of  the  bonding  strap  current  di 
to  an  impulsive  external  field  HQ6(t).  The  normalizing  factor  is  VQ/Rg 
The  numbers  in  parenthesis  indicate  the  values  the  curves  intercept  the 
ordinate  axis. 


TABLE  3.  PEAK  VALUES  FOR  Ife  AND  Ife 


Tl/T2 

0.01 

0.1 

0.2 

0.5 

1 

^5 

T2Ife(peak) 

H  a.sin<fi 
o  2  o 

3.1 

1.57 

1.11 

0.61 

0.36 

0.43  t2 

Tb 

f2ib(peak) 

H  a~sin4 
o  2  o 

370 

37.0 

19.0 

7.4 

3.7 

3.7  t2 

Tb 

(peak)  *  ^  Hosin<j>o 

H^1)  (peak)  »  Hosin4>o  (51) 

T2 


That  part  due  to  direct  field  penetration  can  be  read  off  from  Figure  4  for 
a^l 'a^  *  0.9  and  is  given  by 


u(0),  .  \  „  0.4  „ 

Hi  <Peak>  58  rq  Ho 

•(0) ,  ,  .  3.5  „ 

(peak)  *  — j-  Hq 


(52) 


The  total  peak  penetrant  fields,  (peak)  and  H^(peak),  with  bonding  straps 
are  the  sum  of  Equations  51  and  52.  It  can  be  seen  that  the  bonding  straps 
will  increase  H^peak)  by  a  factor  of  two  and  Hi(peak)  by  as  much  as  an 
order  of  magnitude. 
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VII.  CONCLUSIONS 


It  is  found  that  the  inductive  interactions  among  the  shields  and  the 
presence  of  conducting  straps  bonding  the  shields  reduce  the  effectiveness 
of  a  multisurface  shielded  enclosure  against  the  penetration  of  external 
magnetic  fields . 

For  a  two-surface  spherical  shielded  enclosure  the  important  findings 
can  be  summarized  as  follows.  Let 

p  *  peak  interior  (penetrant)  field 

p  *  peak  time  rate  of  change  of  H^(t)  =  peak  electromotance  force 
(emf)  density 

|HjJ  ■  frequency  spectrum  of  H  (t) 
and  let  the  ratio  of  radii  of  inner  to  outer  shield  *0.9. 

1.  Neglect  of  interaction  of  electrically  thin  shields  underestimates 

•  H,  by  20% 

i.P 

•  H.  by  a  factor  of  4 

i,P 

•  IhJ  by  one  order  of  magnitude  for  f  >  100  kHz  and  enclosures 

with  L/R  time  constant  “  100  us. 

2.  Neglect  of  interaction  of  electrically  thick  shields  under-estimates 

•  |H± |  by  a  factor  of  3  for  f  >  160  kHz  and  enclosures  with  L/R 

time  constant  *  100  us. 

3.  Two  bonding  straps,  each  subtending  a  22.5°  angle  at  the  center 
and  with  inductance  2L^  ■  0.01  L2  (Lj  “  0.2  uH  for  an  enclosure 
of  one-meter  diameter),  increase 

•  H^  p  by  a  factor  of  2 

•  H^  p  by  one  order  of  magnitude 
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APPENDIX  A 
THE  CONSTANTS 


The  constants  A,  B,  C,and  D  that  appear  in  Equation  1  are  given  by 


H  s(t1+ci3t,)  +  s2t,t,(1  -  a3) 

.  o  x _ 4 _ i.  i- _ 

A  ®  y  2  3 

(1+  ST^)  (1+  st2)  -  s  T 2<i 


1  +  ST, 


B  =  -  H 


0  (1  + stl)  (1  + st2)  -  s2T1T2a3 


C  *  “  T 


ST2a 


2  3 

(1  + st1)  (1  + st2)  -  s 


D  *  -  H 


(1+ st^  (1  +  st2)  ~  s  T-^Cl 


where  a  *  a2^ a\’  T1  *  %alCTl^:i/3,  T2  *  ^oa2a2A2^3  ' 

The  constants  A',  B',  C',  and  D'  that  appear  in  Equation  37  are  given  by 


A '  s _ — 

2F 


(KlPI-9K^)sinh  pl  C05h  P2+(K2p2+9K^7)sinh  P2 


l)  (K2p2‘9K^)sinh  P2  cosh  pl"  (KIpl  +  9K^)sinh  P1 


H 

B'  -  .-a 


cosh  P2+(K2P2  +  9i^)3inh',2 


»o(a2 


2F  V  a 


C'  =■-—(—/  (  K,p0  -  _  jsinh  p2 


1/  V  2"2  9K,p 


2P2  / 


where 


s  sp2°2^2 


K 


Uoal 


1  3y1A1  ’ 


K, 


uoa2 


2  3y2A2 


=  cosh  Pl  +(kiPi  +  9K^)sinh  Pl  cosh  p2  +(k2p2  +-^-)sinh  p2 

/a  \3  , 

(9^;  -Kin)(?K57-K2!'2)si,,ht'lSlnh''2 


Unless  and  K2  are  much  smaller  than  unity  one  may  neglect  terms  involving 
the  reciprocal  of  K^p^  and  K2P2  and  obtains  for  F  the  following  accurate 
expression 


F  -  (cosh  p^  +  K^p^sinh  p^) (cosh  p2  +  K2P2sinh  p2>  - 


-  (a2/a1)  K1P1K2p2sinh  p^inh  p2 
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APPENDIX  B 


CALCULATION  OF  '4>Qa£ 

The  magnetic  field  at  the  equatorial  plane  (9  *  tt/2)  can  be  found  from 
Equation  1,  namely. 


H0  =  B  +  C  a3/r3,  a2  <  r  < 
=  D  &2  >  r 


The  magnetic  flux  can  then  be  calculated  via  (Fig.  Bl) 

OAB  =  "  Uo  j  **0 
OAB 

a,  a* 


4> 


-  v 


*?o 


3' 

a. 


^D  +  J  U,+Cl 


rd<j>dr 


-  V, 


2, 
a„<p 
2  o 


'a  '0  \  r 

2 


a.  r 


f  1  -l(  alsin^  \ 

D  +  j  sin  ^  J-r 

a2  L 


B  +  C  -i  jrdr 


-  y. 


2  o 


D  +  B 


ala2 


2  Sin  *o  "T  *o 


+  C  a: 


|>  cosy  -  cos(y  +  $  ) 
o  o 


La2 


a^siny 


1  4-sT 


FH 


(1  +  ST^) (1  + ST2)  -  a3S2T^T2  ^ 


where 


F  -  yoala2S*n<^o 

Tc  ■  t2  J"  (1  -  a)  +  a  tan($Q/2)  (cos<f>0  -  a)/sin$o 
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